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VIBRATION OF A CANTILEVER BEAM THAT 



SLIDES AXIALLY IN A RIGID FRICTIONLESS HOLE 



INTRODUCTION 

Consider a cantilever beam that fits snugly into a frictionless hole, 
Fig. 1. Axes x,y are fixed. A force F(t) moves the beam axially. Hence, 
the length L of the beam outside of the hole is a function of time t. The 
lateral deflection of the beam is y(x,t), where for x < 0, y = 0. The 
axial velocity of the beam is L = dL/dt. 

By Hamilton’s principle, 

(6T + 6W) dt = 0 (1) 

o 




where T denotes the kinetic energy of the system and W denotes the work 
of external and internal forces. The kinetic energy is 



T 




2 , 
dx 



( 2 ) 



where the subscript t denotes partial d i f f eren t i at i on , p = the mass of the 
beam per unit length = constant, and L^ = the length of the beam when t = 0. 
The strain energy of the beam (the negative of the work of internal forces 
for a conservative system) is 



1 



L 




( 3 ) 



where E = the modulus of elasticity of the beam and I = the second moment 
of the beam cross-sectional area. Hence, the virtual work 5W is 

6W = F5L - 5U W 



where 5U = the first variation of the strain energy (Eq. 3)- Because of 
the inertia, the axial movement of the beam tends to bend the beam by beam- 
column action. This effect is neglected here, but it could be included. 
Equations (1), (2), (3), and (k) yield 



[pL^L6L + j py^^(L.t) 5L + F6L - ^ Ely^^(L,t) 6L 



+ P / y^<5yj.dx - El / y^^5y^^dx] dt = 0 



(5) 



If L(t) is considered as given, the terms with the factors 5L and 6L merely 
determine FCt) to give the prescribed motion L(t). The Euler equation for 
the lateral motion of the beam is 



Ely + py 
xxxx tt 



= 0 



( 6 ) 



The natural boundary conditions are 



(7) 
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The forced boundary conditions are 



y(o,t) = y^(o,t) = 0 (3) 

Adopting Hamilton's viewpoint, we may suppose that y(x,o) and y(x,t^ 
are prescribed (Fig. 2). The problem is to determine y(x,t) in region R 
(Fig. 2) subject to the differential equation (Eq. 6) and the indicated 
boundary conditions (Fig. 2). 

The curved boundary at the right side of region R is troublesome. It 
is possible to transform coordinates so that this boundary becomes straight. 

Let 



X = xCc,n) , t = t (c ,n) 

Regard y as a function of (q,n). Then, 

y = y(C,n), C = 5(x,t), n = n(x,t) 

and 



(IV) 
^3x^ t 

‘It) 

d t X 



(^>x 

3t 



{h.1 (ii) + (^) (^) 

3x''t ^3n^c 3x^t 

(1^) (ii) + (ii) (in.) 

^3C^n 3f'x ^3n ? 3t\ 



= y 






3Y\ r3riN 



Here, we regard y as being a function of CC>n) • 
Now let 



c, = , t = n 



(9) 



(10) 



(n) 



0 2 ) 



3 



3 y 3 y / ^ 

When we write , etc., we imply y = y[^,r]) 

and simi larly. 



3y 

Hence, ^ 



(IV) (1) 



Z| 4 

8 y 1 8 y 

— ZT = T * —TT 

8x L 8C 



( 13 ) 



Now, 



and since 



(^) 



= il + ii = 
8C ^at^x 8n 



y 



/ aC '\ _ xL _ CL / a^N _ 
^8t\ L ’ ^at^x 



{]k) 



we find 



and hence, 



y . -(|X)„ f * (fe, 

oC n L art C 



(|f)„ 

8C n 



(^) 

^an''c 



.(^bt) iL - k + 

L ^ac^n L acan 

a% CL ^ + (^) §kl 

acan l ^ac n l ^ac^n l 



(15) 



+ 




Thus, by the third of Eqs. (11), Eqs. (12), Eqs. (14), and Eqs. (15), we 



f i nd 



4>x - 1-4) ^ 

at 8 c 



(ly)k + (- kk) 

^8ca 8can^ ^ L ’ 



r-fi^'iik - (lY^kk + zly.'ikid + 

^ ^8can'^ L ^8C^ L ^8C^ 9n^ 
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and since n = t , we obtain 



2 22*2 '2 2 

fl_y ^ 5_k_ + 9 iil 9 l_y_ ik 

|_2 3C '- 



lY ik 
as L 




in which y = y(^,t) on the right hand side of Eq. (16). 
Consequently, Eqs. (6), (13) and (16) yield 



!■ A + p[i!y i!ki + 2 ^ ^ 



- 9 k_x_ ik . 3y ik ^ £y . ^ 

scat L 3C L g^2^ 



( 16 ) 



(17) 



where y is now regarded as a function of ' and t. Equations (7) and (8) 
become 



yrr (' >t) = .t) = 0 

(13) 

y (0,t) = y. (0,t) = 0 
>- 



Equation (17) is much more complicated than Eq . (6). However, the 
region R is simpler (Fig. 3)- The problem remains a linear boundary-value 
problem if in Eq. (17) we regard L(t) as a given function of time t. If 
we consider L(t) as unknown, we must obtain another equation which defines 
L(t). This equation is obtained from Eq . (5) from the terms with the 
factors 6L and 6L. Hence, we find 



5 






( 19 ) 



In terms of (C,t). Eq- (19) becomes 



L + Yrcd .t)- p(“ f Yrd .t)+ y^d.t)]^) 

p o L s 



(20) 






F(t) 



Equation (20) is subject to the initial conditions 

L(0) = 

L(0) = 0 



( 21 ) 



where is the initial length of the beam at time t = 0, (see Fig. 1). 

Thus, if y(x,t) and L(t) are considered as the unknowns, Eqs. (17) 
and (20). subject to the boundary conditions of Eqs. ( 1 8) and (21), respec- 
tively, represent a complicated nonlinear boundary value problem. 

ANALYSIS OF EQUATIONS 

(a) L(.t) Prescribed. If L(t) is given, then the problem reduces to 
seeking a solution of Eq. (17) subject to the boundary conditions, Eqs. (18). 
The problem is still of considerable difficulty. However, the problem is 
a linear boundary-value problem. The nature of the solutions of Eq . (17) 
is unknown. indeed, one must consider whether or not a unique solution of 
the boundary-value problem exists. Also, is there a unique solution of 
the correspond i ng initial value problem in which the condition y = F^ (^) 
on t = t^ is discarded and the functions y(5,0) and y^(5,0) are prescribed? 
One way to approach such questions is through the calculus of finite 



6 



differences, in which the differential equation is replaced by a linear 
partial difference equation. The existence of a solution of the difference 
equation is a question of linear algebra. The most practical way of 
developing numerical solutions of Eqs. (17) and (18) is an open question. 

One possible approach is through the calculus of finite differences, or a 
combination of the calculus of finite differences with the Runge-Kutta 
method. Another approach is through finite-element methods. 

The problem is simplified greatly if the axial motion is so slow 
that time derivatives of L(t) are negligible. For example, if the beam is 
vibrating in a natural mode and it is then drawn slowly into the hole, 
does it continue to vibrate in a natural mode? How do the frequency and 
amplitude vary with L under these conditions? For a rigid, frictionless 
hole, the total mechanical energy is constant, since no energy is imparted 
at the root. This criterion serves to determine how the amplitude varies. 
However, we should not have to apply this criterion, since constancy of 
mechanical energy should follow automatically from the solution of Eq . (17)* 
[f L and L are negligible, Eq. Cl7) reduces to 






( 22 ) 



Equation (22) is much simpler than Eq. (17)* However, the fact that L is 
a prescribed function of t still complicates it. There is a question as 
to whether the character of the solution is altered by the reduction of 
Eq. Cl7) to Eq. (22). Again, does the initial value problem posed by 
Eqs. (18) and (22) have a unique solution? If so, the constancy of mechanical 
energy, as a separate principle, is not needed to determine the dependency 
of amplitude upon L. In fact, constancy of mechanical energy should be a 
deducible consequence of Eqs. (18) and (22). 
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Unlike Eq . (17), Eq. (22) allows a theory of natural modes. For 
example, consider solutions of Eq. (22) of the form y = f(q)g(t). Then, 



Eq. (22) yields 



Hence, 



f E I g 



f" _ pl'* g" _ .k 

f ^ ’ eT i ^ 



in which g is a constant. Then, 



g 



1 1 




= 0 



Unfortunately, Eq. (24) does not yield simple harmonic motion, 
is a function of t. However, the equation 



II 1 




= 0 



can be integrated. With the root conditions (Eqs. 18) f(o) = f 
Eq. (25) yields 

f(C) = A(sinhBC * sinBC) + B(cosh6C - cosBC) 

Now, the free-end conditions (Eqs. 18) f"(l) = f'"'(l) = 0 yield 

A(sinhS + singL + B(coshB + cosB)_ = 0 
A(cosh6 + cos6)_ + B(sinhB - sinB) = 0 

and the vanishing of the determinant of Eqs. (27) yields 



(23) 

(24) 

because L 

(25) 

’ (o) =0, 

( 26 ) 

(27) 
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coshB cosg = -1 



(28) 



These results are exactly like that for natural modes of a cantilever with 
fixed length. The roots of Eq. (28) are constants 6 | , 62 * 6 ^. ••• • They 
have been tabulated by Young and Felgar (1). 

As is customary, let B = 1 and A = -a^ . Then, by Eq . (27), 

coshB + cosB 

a = ^ ( 29 ) 

n sinhB + sinB 
n n 

Consequently, , ••• are constants, independent of L. The values of 

have been tabulated by Young and Felgar. Equation (26) gives the natural 
modes , 



f (g) = coshg g - cos6 ^ - a (sinhg g - s i n6 C) (30) 

n n n n n n 

The functions ^ ^fe orthogonal in the interval (0,1); in fact, it is 
eas i 1 y shown that 



1 




o 



f f dC 
n m 




f "f "dC = 0 , m n 
n m 



1 




o 




(f '')^dc = 

n n 



(31) 



The fact that the functions f^(C) satisfy the end conditions (Eq. l 8 ), 
makes then ideal approximating functions for use in numerical attacks on 
Eqs. (17) and (22). Presumably, a solution of Eqs. (17) or (22) can be 
approximated by 
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( 32 ) 



N 

y(5,t).^ q„(t)f„(S) 

n=l 



This approach reduces the problem to one of a finite number of degrees of 

freedom. The generalized coordinates are q (t). The Lagrange theory of 

n 

linear vibrations is then applicable. Accordingly, let 



03 



n 




( 33 ) 



Then, Eqs. ilk) becomes 

g ’• + 03 ^g = 0 
n n ^n 



( 34 ) 



Since L depends on t, note that 03^ depends on t. Consequently, 03^ is not 
exactly the frequency of a harmonic motion. The general solution of Eq. 
(34) is 

9 „ - < 35 ) 



in which A , B are constants of integration. An infinite series solution 
n n ^ 

of Eq. (22) is then 



y(5.t) 

n = l 



f iOgit, A B^) 

n n n n 



( 36 ) 



Equation (36) au toma t i ca 1 1 y satisfies boundary conditions (Eqs. 18) at the 
free end and the root. One might expect that the constants (A^ , B^) can be 
chosen to satisfy the conditions at t = 0 and t = t^; i.e., 



1 0 



CO 



L ‘n- “n> ‘ 

n=l 



(37) 



2] f„‘5>9n"r \’®n> - 

n=l 



Alternatively, the constants A , B might be chosen so that the beam has 
a given initial deformation (the first of Eqs. 37) and given initial 
velocities; i.e., y^(C»0) might be specified instead of y(C,tp. 

If L(t) is a sufficiently slowly varying function of t, the approxi- 
mate solution of Eq. (3^) is 



g = A sinoj t + B cosoo t 
n n n n n 



Then, Eq. (361 becomes 



y(C.t) = 






(A. 



s i noj t + 
n 



B cosoi) t) 
n n 






n=l 



The initial value problem then reduces to 



(38) 



(39) 



V B f (C) = FiO 
/ ^ n n o 



n=l 



(40) 



Z A w °f (0 = 4>,(a 

n n n o 



n=l 



in which d) (E) = y (q,0), and OJ is the initial value of (jO . Under fairly 
o t n n 

broad conditions, the constants A , B can be chosen to satisfy Eq. {kO) 

n ’ n 

when F’q(^) and <^q('^) are prescribed functions. Then Eq. (39) represents the 
solution of the initial value problem. 



1 1 



The kinetic energy and the strain energy are, respectively, 



I 



T = 




00 







2 2 
oj (A COS03 t - B sino3 t) 
n n n n n 



1 



U 




2 E 1 \" > L 

y--d5 = — ^ / B (a sin03 t + B cosoo t) 

CC 2L^ n n n n n 



Hence, in view of Eq. (33), the total mechanical energy is 



T + u = ^ 6 ^(A ^ + B 

2|_3 n ' n n 

n=l 



Thus, the total mechanical energy is independent of t, provided that A^ 

3/2 

and B^ vary as L . Apparently, for the conditions stated, this require- 
ment ensures conservation of mechanical energy. 

(b) F(t) Prescribed. If one prescribes F(t), then Eqs. (17) and (20) 
subject to Eqs. (18) and (21) must be solved simultaneously in order to 
determine L(t). The method of solution for both cases L(t) prescribed 
and F(t) prescribed is outlined in the following section. 



Method of Solution 

Equations (17) and (20), together with the initial and boundary condi- 
tions given by Eqs. (18) and (21) form a system of two nonlinear partial 
differential equations in the two unknowns L(t) and y(x,t). An approximate 
numerical solution of these equations can be obtained by a Galerkin finite 



e 1 emen t f o r mu 1 a t i on . 



Thus, in accordance with the Galerkin FEM, we form the approximate 



solutions of L and y as (\> and v which are, 



L(t) = (})(t) = 6(t) 
y(x,t) = v(x,t) = N^(x)a(t) 



(41) 



where N are a set of quadratic shape functions with local support. A 
quadratic shape function is associated with each nodal point of the dis- 
cretized finite element model. Quadratic Interpolation over an element 
requires each element has three nodal points, say one at each end of the 
element, and a nodal point in the center of each element. 

Figure k shows that the odd numbered shape functions , ..., 

associated with nodal points at the ends of an element span five nodal 
points (two elements), while the even numbered shape functions ..., 

span the three nodal points of one element. 

In accordance with the Galerkin FEM, we form the residual function 
for Eqs. (I7l and {10) as. 



* 2 -2 

R (x.t) = ^ v"“ + ^ + 2v'r ±- 

<p <p 0 



2v'C I - v'5 I + v] 

V P 



(42) 



1^2 ~ 2qL ~ p{--^v'(1,t) +v(l,t)} ] 

(j) ^ 






F(t) 



(43) 



The finite element equations are obtained by requiring that the residual 
functions, given by Eqs. (42) and (43), be orthogonal to each of the basis 
functions. That is. 
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1 

J NR^dx = 0 (44) 

0 

1 

y* HR^dx = 0 (45) 



where H is the function given by 



H = 1 0 < X < 1 



(46) 



Substitution of Eq. (41) into (42) and (43) gives 



• 2 * 2 

R (x,t) = -!j-(E I N"''’ a) " + ^ N"''’a + 2? N'''’a 

B B B 



(47) 



- 2 ? I N'"'’a - 5 -I N'''’a + N''’a] 



o 6 



C=1 



p{- I N'"''a 



+ N^d 



?=1 






C=1 



o 



(48) 



Substitution of Eqs. (47) and (48) into Eqs. (44) and (45) respectively, 
and performing integration by parts, gives 



B 



2 T 

M » M I ' r 



B 



BT + / EIN‘‘N" a^a- ^ N ’ N ' 0§a + 2 / ^NH^ d^a 






e 



e 



6 



J.r., _ 3 



2 -| / CNN' dCa--|; ^NN' dCa + / NN d?d = 0 



(49) 
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and 



S + ■TT'i 7T [ ( T , ~ + :r , ) 

2pLq ^ ^2 2n-l ^ 2 2n ^2 2n-l 

e e e 



( 50 ) 



/ B/ 1 
P{" a 



^2n-l ■ 2'^ “2n ^2n+l ) + ^ 



where BT denotes a boundary term. 
Lett i ng , 



1 




e e e 



( 51 ) 
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Equations (^9) and (51) become 



, S 



• 2 *2 

BT + Aa - Aa + 2 ^ Ca - 2 Cot + Da = 0 



3 + 



B 

El 



^^[3-pb]=^P(t) 
o B o 



( 52 ) 



(53) 



Defining the following terms, 
Y = -(BT + Aa) 



0 ) = 



C = 



1 • ' ' R 

- — ^ B Aa + 2 — j Ca 

B ~ B 



B 



-2 C 

B ^ 



(5A) 



! = ■ B 



G = F(t) - (a - pb) 



pL ''' 2pL 
o o 



where B is the value of B at the previous time, Eqs. (52) and (53) become. 



Da + eB + C a + o)B 

B 

or, in matrix form 

1 2n+l 2n+2 



2n+l 

2n+2 



Y 

G 



0 . 



(55) 




solved by numerical integration. 
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Figure 2. Region of Integration for Equation (6). 
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Figure 3- Region of Integration for Equation (17) 



* 6i.eM.£N>p pqiHT 
XfcL£M. P0»HT 



£NP ■ gUM. MO^AW PqiWT 
GNP-GUOBAU NOSAU PoiMT 




Figure h. Discretized FEM Model with Quadratic Basis Functions 
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